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ABSTRACT 

Context. Large-scale magnetic fields resulting from hydromagnetic dynamo action may differ substantially in their time dependence. 
Cyclic field variations, characteristic for the solar magnetic field, are often explained by an important fi-effect, i.e. by the stretching 
of field lines due to strong differential rotation. 

Aims. The dynamo mechanism of a convective, oscillatory dynamo model is investigated. 

Methods. We solve the MHD-equations for a conducting Boussinesq fluid in a rotating spherical shell. For a resulting oscillatory 
model, dynamo coefficients have been computed with the help of the so-called test-field method. Subsequently, these coefficients have 
been used in a mean-field calculation in order to explore the underlying dynamo mechanism. 

Results. The oscillatory dynamo model under consideration is of a 2 £l-type. Although the rather strong differential rotation present in 
this model influences the magnetic field, the fl-effect alone is not responsible for its cyclic time variation. If the £l-effect is suppressed 
the resulting a 2 -dynamo remains oscillatory. Surprisingly, the corresponding trfi-dynamo leads to a non-oscillatory magnetic field. 
Conclusions. The assumption of an fffl-mechanism does not explain the occurrence of magnetic cycles satisfactorily. 



1. Introduction 

The study of the solar cycle has motivated dynamo theory for 
many decades. Hence, the solar dynamo has become the pro- 
totype of oscill atory dynamos. H owever, its explanation is still 
controversial dJones et alj l2010h . Most solar dynamo models 
have been built on the a ssumption of an ffQ-dynamo mecha- 
nism (Ossen driiverll2003h : that is, the poloidal field results from 
the interaction of helical turbulence with the toroidal field (a- 
effect) whereas the toroidal field is thought to originate from the 
shearing of poloidal field lines by strong differential rotation (fl- 
effect). This attempt is attractive for mainly two reasons: 

First, the existence of a strong shear layer at the bottom of 
the solar convection zone is observationally well established, and 
the importance of a resulting Q-effect i s non-controv ersial. 

Second, Parker's plane layer mo del (|Parkerll95 5 |) and i n par- 
ticular mean-field electrodynamics (Steenb eck et al.l ll966) pro- 
vide a very elegant theoretical framework for this approach. 
Within mean-field theory, attention is focused on large scale, 
i.e. averaged fields, only, and the induction equation may be 
replac ed by a mean-field dynamo equation (iRrause & Radlerl 
1980) 



dB - - - 

— = Vx(£+VxB-nVxB), 
at 



(1) 



in which B and V denote the average magnetic and the average 
velocity field, rj stands for the magnetic diffusivity and £ is the 
mean electromotive force. Moreover, it is assumed that £ is ho- 
mogeneous in the mean magnetic field and may be replaced by 
a parameterisation in terms of B and its first derivatives 



£ = aB + bVB . 



(2) 



In 0, the so-called dynamo coefficients a and b are tensors of 
second and third rank, respectively, and depend only on the ve- 
locity field and the magnetic diffusivity. The traditional a-effect 



implemented in a large number of solar dynamo models (e.g . 
[ Steenbeck & Krausdll969l:lRobertslll97llRoberts & Stixlll972b 
Stix 1976; Ossendriiver 2003; Brandenburg & Subramanian 
2005; IChan et aTl |2008) corresponds to the isotropic compo- 
nent of a in relation (0, while the Q-effect results from the <p- 
component of the Vx(FxB) term in equation (0. 

However, a strong differential rotation is not a neces- 
sary condition for oscillatory solutions of the dynamo equa- 
tion (0, as has been demonstrated in several papers (see e.g. 
Ritdler & Brau erlll987l ISchubert & Zhangl 12000: iRiidiger et all 
2003; Ste fani & Ger beth 2003). These authors construct models 
in which the toroidal field is likewise generated from the poloidal 
field by an a-effect (a 2 -models) and investigate necessary con- 
straints on a, the boundary conditions for the magnetic field and 
the geometry of the dynamo region in order to obtain oscilla- 
tory solutions of 0. Recently, oscillatory dynamo models have 
also been in vestigated by means of direct numerical simulations. 
Mifr a et a"D (1201 Oh performed dynamo simulations in a wedge- 
shaped spherical shell with an applied forcing and demonstrate 
again the existence of oscillatory or 2 -dynamo models. 

However, the success of mean-field models in reproducing 
solar-like variations of the magnetic field relies partly on the 
large number of free parameters, i.e. on the arbitrary determina- 
tion of the dyna mo coefficients a and b. An alternative approach 
is presented by iPetrelis et al.l d2009l) . They construct amplitude 
equations guided from symmetry considerations and analyse po- 
larity reversals and oscillatons of the magnetic field resulting 
from the interaction between two dynamo modes. 

Self-consistent, global, convective dynamo models with 
cyclic magnetic field variation s have been bubl is hed b y 
iBusse & Simitevl (120061) . and iGoudard & Dormvl J2008). 
Convective dynamo simulations with stress -free mechanical 
boundary conditions dBusse & Simitevl 120061) exhibit a strong 
and a weak field branch, depending on the initial conditions for 



2 



M. Schrinner et al.: Oscillatory dynamos and their induction mechanisms 



the magnetic field. If the magnetic field is initially weak, stress 
free boundary conditions enable the development of a strong 
zonal flow carrying most of the kinetic energy and rendering 
convection ineffective. The magnetic field resulting from these 
dynamos is rather small scaled, often of quadrupolar symmetry 
and weak. Oscillatory solutions of the induction equation are 
typical for this dynamo branch. 

A transition from steady to oscillatory dynamos may 
also be governed by t he width of the convection zone; 
iGoudard & Dormvl d2008f) found oscillatory models by decreas- 
ing the shell width. In this study, we follow their approach and 
analyse the dynamo mechanism for these oscillatory models. In 
particular, we address the question whether an Q-effect is re- 
sponsible for the cyclic variation of the magnetic field. Different 
from previous work, we determine the dynamo coefficients a and 
b from di rect numerical simulation s with the help of the test-field 
method (ISchrinner et al.l2005ll2007l) . The application of a and b 
in a mean-field calculation reveals their importance for the gen- 
eration of the magnetic field. 

2. Dynamo calculations 

We consider a conducting Boussinesq fluid in a rotating spheri- 
cal shell and solve the equations of magnetohydrodynamics for 
the velocity v, magnetic fi eld B and temperature T as given by 
IGoudard & Dormvl d2008l) with the help of the code PaRoDy 



in which the linear operator D is defined as 



DB = V x B + aB 



bVB - —V x B . 

Pm 



(7) 



The time evolution of each mode is determined by its eigenvalue 
cr and proportional to exp(crf). Fo r more details concernin g the 
eigenvalue calculation, we refer to lSchrinner et al.l d201 0b). 

We also consider the evolution of a kinematically advanced 
magnetic field, Bj r , governed by a second induction equation 



dB Tr 1 , 

— — ^ = V x (v X Bj r ) + — — V Bxr • 
Of Pm 



(8) 



The tracer field Bj r experiences the self-consistent velocity field 
at each time step but does not contri bute to the Lorentz force 
and is passive in this sense (see also lSchrinner et aI1l2010ah . Its 
evolution will be compared with mean-field results originating 
likewise from a kinematic approach. Moreover, a kinematically 
advanced tracer field allows us to test for the influence of the 
Q-effect in direct numerical simulations. In a numerical experi- 
ment, we subtract the contribution of the O-effect and the mean 
meridional flow in the equation for the tracer field, 



dB 1 



1 



dDormv et al.l dl998l) and further developments), 



V x (v x B Tr ) + — — V Brt - V x (V x B Jr ) , (9) 
of Pm 

and study in this way the outcome of a kinematic o^-dynamo. 



E ' + V VV ~ V2v ) + 2Z X V + V ^ = 

r 1 
Ra — T + (Vxfi)xB, 

r„ Pm 



dB 



1 



— = V x (v x B) + V Z B , 

dt Pm 

5 + (vV)(r + r s ) = iv 2 r. 

of Pr 



(3) 
(4) 
(5) 



Governing parameters are the Ekman number E = v/Q.L 2 , the 
(modified) Rayleigh number Ra = argoATL/vQ., the Prandtl 
number Pr = v/k and the magnetic Prandtl number Pm = vjr\. 
In these expressions, v denotes the kinematic viscosity, Q. the 
rotation rate, L the shell width, aj the thermal expansion coef- 
ficient, go is the gravitational acceleration at the outer boundary, 
Ar stands for the temperature difference between the spherical 
boundaries, k is the thermal and 77 = 1/ficr the magnetic diffusiv- 
ity with the magnetic permeability fi and the electrical conduc- 
tivity cr. Furthermore, the aspect ratio is defined as the ratio of 
the inner to the outer shell radius, r,/r ; it determines the shell 
width. 

In our models, convection is driven by an imposed temper- 
ature gradient between the inner and the outer shell boundary. 
The mechanical boundary conditions are no slip at the inner and 
stress free at the outer boundary. Moreover, the magnetic field is 
assumed to continue as a potential field outside the fluid shell. 

Time-averaged dynamo coefficients for an axisymmetric 
mean magnetic field have been determi ned from direct numer- 
ical simulations as described in detail by ISchrinner et alj d2007l) 
and as rec e ntly discussed for time-dependent dynamo models by 
Schrinner (2011). In a second step, these coefficients have been 
applied in a mean-field model based on equation ([TJ written as 
an eigenvalue problem, 



3. Results 

The model under c o nsid eration has been previously studied by 
IGoudard & Dormvl (120081) . It is defined by E = 10 3 , Ra = 
100 (= 2.8 Ra c ), Pm = 5, Pr = 1 and an aspect ratio of 0.65. 
Except for the stress-free mechanical boundary condition ap- 
plied at r — r a and an increased aspect ratio, the governing pa- 
rameters are those of a rather si mple, quas i-steady benchmark 
dynam o (jChristensen et ai] |2001l) . However, Goudard & Dormy 
(2008) report a transition from steady, dipolar to oscillatory 
models for these parameter values. Note that the model re- 
quires a rather high angular resolution up to harmonic degree 
/ = 1 12 

(■max — 11^.. 

Figure Q] displays the radial component of the velocity field 
at a given radial level. A typical columnar convection pattern 
is visible, even though the convection columns are noticeably 
disturbed by the influence of the curved boundaries and a strong 
zonal flow carrying about 50 per cent of the kinetic energy. The 
magnetic Reynolds number based on the rms-velocity and the 
shell width, Rm = v^L/ri, is about 90. The flow is symmetric 
with respect to the equatorial plane and convection takes place 
only outside the inner core tangent cylinder. 

The evolution of the magnetic field is cyclic. In figure|2](top), 
contours of the azimuthally averaged radial magnetic field at 
the outer shell boundary varying with time are plotted in a so- 
called butterfly diagram. A dynamo wave migrates away from 
the equator until it reaches mid-latitudes where the inner core 
tangent cylinder intersects the outer shell boundary. The mag- 
netic field looks rather small scaled and multipolar. This is con- 
firmed by the magnetic energy spectrum which is essentially 
white, except for a negligible dipole contribution. Furthermore, 
the magnetic field is weak, as expressed by an Elsasser number 
ofA = 5^/0^^ = 0.13. 

The kinematically advanced tracer field grows slowly in 
time, i.e. the model under consideration is kinematically unsta 



crB = V x DB , 



tcai iy 1 

(6) ble according to the classification by ISchrinner etal.1 feoiOal) . 
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Fig. 1. Snapshot of the radial velocity of the considered dynamo 
model at r — 0.79 r . The velocity component has been nor- 
malised by its maximum absolute value, v rimax = 24.46 v/L. 
Hence, the colour coding ranges from — 1, white, to +1, black. 
Contour lines correspond to +0.2 and +0.6. 



But, deviations of the tracer field from the actual field are hardly 
noticeable in the field morphology. Moreover, the very same 
dynamo wave persists in the kinematic calculation (see also 
iGoudard & Dormvl d2008l) ). as visible in figure [2] (middle). Note 
that the tracer field in figure [2] has evolved from random initial 
conditions. 

A mean-field calculation based on the dynamo coefficients 
a, b and the mean flow V determined from the self-consistent 
model is presented in the bottom line of figure [2] The fastest 
growing eigenmodes form a conjugate complex pair and give 
rise to a dynamo wave which compares nicely with the direct 
numerical simulations. Since this model depends on the full a- 
tensor and the mean flow, we refer to it as an a 2 Q-dynamo. 

The influence of the differential rotation may be suppressed 
in the kinematic calculation of the tracer field without changing 
any other component of the flow. A butterfly diagram resulting 
from a kinematically advanced field according to equation (0 is 
presented in figure [3] (top). The evolution of the magnetic field is 
again cyclic. Apart from small-scale variations on shorter time 
scales, a dynamo wave migrates from mid-latitudes towards the 
equator. This is in agreement with a corresponding mean-field 
calculation in which the mean flow V in (O has been canceled: 
The bottom chart of figure[3]provides the butterfly diagram stem- 
ming from the fastest growing eigenmodes of the resulting a 2 - 
dynamo. An explanation why direct numerical simulations and 
mean field calculations compare somewhat better in figure|2]than 
in figure[3]is provided in appendixlAl 

The time evolution of the related aQ-dynamo is of further in- 
terest. As the a-effect is not directly accessible in direct numer- 
ical simulations, the corresponding ffQ-dynamo can be realised 
in a mean-field calculation, only. In a first attempt, we have set 
a rr = &ee = in order to suppress the generation of toroidal 
field from poloidal field by an a-effect. Both components make 
major contributions to this process. The leading eigenmode re- 
sulting from this calculation is shown in figure |4j it is real, i.e. 
non-oscillatory, and close to marginal stability. The results re- 
main similar if we neglect further, non-diagonal components of 
a. 

4. Discussion 

The Frequency and the propagation direction of the dynamo 
wave visible in figure [2] d epend strongly on the di fferential ro- 
tation, in agreement with iBusse & Simitevl (120061) . We follow 
their approach and give an estimate fo r the cycle f requency ap- 
plying Parker's plane layer formalism dParkerll955l) . To this end, 
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Fig. 2. Azimuthally averaged radial magnetic field at the outer 
shell boundary varying with time (butterfly diagram) resulting 
from a self-consistent calculation (top), kinematic calculation 
according to ([8]) (middle) and mean-field calculation (bottom). 
The contour plots have been normalised by their maximum abso- 
lute value at each time step considered. The colour coding ranges 
from -1, white, to +1, black. 

we introduce a Cartesian coordinate system (x, y, z) correspond- 
ing to the {(/>, 6, r) directions and define mean quantities to be x- 
independent. Moreover, we write B = Be x + B p — Be x + V x Ae x 
and reformulate (Q]) in the following simplified manner 



dA 
dt 
8B 
dt 



a xx B + — V 2 A , 
Pm 

d dA dV x 8A 1 , 

-7r(ok-7r-) + — - + — v B ■ 
oy oy az dy Pm 



(10) 



(11) 



In the above equations, we have considered only the dominant 
diagonal components of a, a„ and a^ corresponding to a zz and 
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Fig. 4. The leading dipolar eigenmode resulting from a mean- 
field calculation with a rr = agg - 0. Contour plots of all three 
components are presented, each normalised separately by their 
maximum absolute values. Maxima and minima are written next 
to each plot. The colour coding ranges from — 1, white, to +1, 
black, and contour lines correspond to +0.1, +0.3, ±0.5, ±0.7 
and ±0.9. 



U O.Z O.A O.fi £}.B t.q 1.2 
Magnetic dinuston lime 

Fig. 3. Azimuthally averaged radial magnetic field at the outer 
shell boundary varying with time (butterfly diagram) resulting 
from a kinematic calculation with subtracted O-effect (top) and 
a corresponding mean-field calculation (bottom). The contour 
plots are presented as in figure ©. 



a xx ; all components of b and the mean meridional flow have been 
neglected. Furthermore, V has been assumed to depend only on 
z. Then, the ansatz 



(A, B) = (A, B) exp(/jfc • x + erf) , 



leads to 

pA = a xx B, 



P B = 



a zz k 2 + iky 



( dV da z: 
dz dy t 



A, 



(12) 

(13) 
(14) 



with p — cr + \k\ 2 /Pm. From dl3K JT4l . we derive a dispersion 
relation 



dz dy t 



(15) 



from which the real and the imaginary part of cr can be calcu- 
lated. If a xx is positive (e.g. in the northern hemisphere), it fol- 
lows 



A = %(cr) = -\k\ 2 IPm 



y (a zz k 2 ,) 2 + 



'dV da zz 
, dz dy , 



kl + a 7 ?k?. 

y iz y 



1/2 



(16) 




— \ (a zz k 2 ) 2 + 



'dV da zz 
^dz dy , 



k\ cx 77 k v 



1/2 



. (17) 



The sign in ( ITTb is determined by the sign of k v (dV /dz-da zz /dy). 
If we further assume that the frequency is dominated by differen- 
ti al rotation and neglect th e a-terms in (JTTJ, we estimate similar 
to lBusse & Simitevl (l2006l) 



3(cr) * ± \^a xx ^ T 



1/2 



(18) 



In ( TT8l . Ej denotes the the kinetic energy density due to the ax- 
isymmetric toroidal velocity field and k v m In/L has been used. 
Approximating a xx by the rms-value of a^, a xx = 11.50v/L, 
and with Et = 71 v 2 /L 2 , we find u> « ±103.7 rj/L 2 which is 
surprisingly close to a> — +lOQ.95rj/L 2 in the full calculation 
presented in figure [2] Note that a xx and dV/dz are of the same 
order of magnitude and contribute equally to u>. The sign in ( TT~8b 
is determined by the sign of the product a^dV^/dr which is 
positive in the northern and negative in the southern hemisphere. 
Therefore, our estimate in ( TT8b predicts a dynamo wave migrat- 
ing away from the equator. This is in agreement with the simu- 
lations shown in figure [2] 

However, the attempt to describe the model under consid- 
eration as an aQ-dynamo fails. An oscillatory mode with a fre- 
quency close to the above estimate turns out to be clearly subcrit- 
ical in a mean-field calcuation, if a rr and agg are omitted. Instead, 
this model is governed by a real, dipolar mode close to marginal 
stability (see figure|4|. Hence, the Q-effect is only partly respon- 
sible for the generation of the mean azimuthal field, as confirmed 
by figure [5] The chart in the middle compares the Q-effect, 
rB, d(r~ 1 V#) /dr + r~ 1 sin Q Bg d(sin 6~ l V^/dO in grey scale with 
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the mean azimuthal field displayed by superimposed contour 
lines. In particular, the elongated flux patches close to the inner 
core tangent cylinder are, if at all, negatively correlated with the 
Q-effect. Consistent with this finding, the poloidal axisymmetric 
magnetic energy density exceeds the toroidal one by 20%. 

Differential rotation alone is not responsible for the cyclic 
time evolution of the magnetic field, despite its influence on the 
frequency and the propagation direction of the dynamo wave. 
This is most clearly visible in figure [3] Simulations without dif- 
ferential rotation still lead to a dynamo wave even though its fre- 
quency and propagation direction have changed. In the frame- 
work of Parker's plane layer formalism, the frequency of this 
oscillatory <? 2 -dynamo depends crucially on -da zz /dy instead of 
dVx/dz- Note the additional minus sign, which might explain 
the reversed propagation direction if the assume that da zz /dy is 
predominantly positive. But different from the radial derivative 
of the mean azimuthal flow, {l/r)da rr /d6 is highly structured, 
changes sign in radial direction and exhibits localised patches of 
low negative values (see figure[5]i. Therefore, we do not attempt 
to give an estimate for the frequency similar to ( TT~8T >. 

In order to better understand the influence of the mean 
flow on the frequency of the dynamo wave, we have gradu- 
ally changed the amplitude of V in a series of kinematic cal- 
culations. Results are presented in figure [6] Stars denote fre- 
quencies obtained from eigenvalue calculations according to (O, 
wheras triangles stand for frequencies estimated from kinematic 
results due to equation (0. In both cases, the amplitude of the 
mean flow V has been varied by multiplication with a scale fac- 
tor /. For / = 1, the original calculation is retained, while for 
/ = 0, we reproduce the a 2 -dynamo already discussed above. 
Frequencies of dynamo waves resulting from direct numerical 
simulations according to © have been meassured for / = 1 , 0.7, 
and 0.5. Owing to the turbulence present in the simulations, these 
are rather rough estimates and error bars have been included. 
Nevertheless, the results obtained are in satisfactory agreement 
with the eigenvalue calculations. The frequencies in figure [6] 
decrease continously with decreasing scale factors. If the am- 
plitude of V is reduced to 25 per cent of its original value, oj 
changes sign and the propagation direction of the dynamo wave 
is reversed. The dashed-dotted line in figure|6]gives o> according 
to (fT8l ) as predicted for an <mi-dynamo. It matches the numer- 
ical results if dV/dz dominates in (fTTT i but deviates clearly for 
smaller amplitudes. On the other hand, it is illustrative to use 
relation ( fTTT i to model the dependence of u on the mean flow. 
If we set da zz /dy - 0.25 dV/dz and determine a representative 
value for a zz inverting ([T71 > for dV/dz = and u = -29.15j]/L 2 , 
the dashed line in figure[6]results from dTvT >- It fits the numerical 
data rather well and converges towards the frequencies predicted 
for an <m>-dynamo, if the amplitude of V is sufficiently high. 

Let us stress again that some caution is needed in applying 
the present mean-field analysis to non-linear direct numerical 
simulations, as the the dynamo model considered here is kine- 
matically unstable. Strictly speaking, our mean-field results are 
only relevant for the kinematically advanced tracer field. But, 
because the model is close to dynamo onset and only weakly 
non-linear, we believe that our interpretation is also valid for the 
fully self-consistent field. This is in particular confirmed by the 
rather good agreement of the three butterfly diagrams presented 
in figure [2] 



dVJdr Q-effect (Mr) dajdd 

Max: 59.46 Max: 1.19 Max: 36.56 

Min : -7.85 Min : -1.19 Min : -49.03 




Fig. 5. Left: dV^/dr in units of v/L 2 . Middle: Q-effect as given 

by rB r d(r l V,p)ldr + r l s,\n6 B B dinner 1 V,p) 1 86 (greyscale) 

and B,f, (superimposed contour lines, solid [dashed] lines indi- 
cate positive [negative] values). Right: (l/r)da rr /d9 in units of 
v/L 2 . The contour plots are presented in the same style as in fig- 
ure g] 



150 ■ 




0.2 0.4 0.6 0.8 1 
scale factor f 



Fig. 6. Frequencies resulting from kinematic calculations in 
which the amplitude of the mean flow has been changed by mul- 
tiplication with a scale factor, /. Stars denote frequencies stem- 
ming from an eigenvalue calculation according to ©, whereas 
triangles are estimates obtained from kinematic results due to 
(O. The dashed-dotted line gives frequencies as predicted for 
an <m>-dynamo by ([T8l >, while the dashed line represents co as a 
function of V for an a 2 w-dynamo according to ( fTTI ). 



5. Conclusions 

A particular dynamo mechanism does not seem to be responsi- 
ble for the occurrence of periodically time-dependent magnetic 
fields. It turns out, that the influence of the large-scale radial 
shear (the Q-effect), is not necessary for cyclic field variations. 
Instead, the action of small-scale convection, represented by a 
spatially structured dynamo coefficient a rr , happens to be essen- 
tial. For the model presented here, small convective length scales 
are forced by a thin convection zone. Further investigations are 
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needed to assess whether our finding is representative for a wider 
class of oscillatory models. 
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Appendix A: The use of time averaged dynamo 
coefficients 

In the following, azimuthal averages are, as throughout in the 
paper, denoted by an overbar, time averages are expressed by 
brackets, < ■ ■ • >. Initially, dynamo coefficients have been de- 
termined for an azimuthally averaged, mean magnetic field B. 
Hence, the evolution of the latter is given by 

c\n r 

— = Vx(aB + bVB + VxB VxB). (A.l) 

ot Pm 

But, the dynamo coefficients a, b and the mean flow V vary 
stochastically in time. In order to describe the average dynamo 
action, we take in addition the time average of these quantities 
and write approximatively, 

8B - - - - 1 - 

— » Vx(<a> B+ <b> VB+ <V> xB VxB). (A.2) 

ot Pm 

We emphasise that there is no a priori relation between the left 
hand side and the right hand side of equation (IA.21 >. The actual, 
azimuthally averaged magnetic field will deviate from our mean- 
field description the stronger, the more a, b and V fluctuate in 
time. Among these three quantities, the mean flow V is almost 
time independent, whereas a and b vary considerably. This is 
the reason, why the butterfly diagrams in figure [2] are in better 
agreement than in figure [3] for which the stabilizing influence of 
the mean flow has been omitted. 
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